Isospin splitting of the nucleon mean field 
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The isospin splitting of the nucleon mean field is derived from the Brueckner theory extended to 
asymmetric nuclear matter. The Argonne V18 has been adopted as bare interaction in combina- 
tion with a microscopic three body force. The isospin splitting of the effective mass is determined 
from the Brueckner-Hartree-Fock self-energy: It is linear acording to the Lane ansatz and such that 
> nip for neutron-rich matter. The symmetry potential is also determined and a comparison 
is made with the predictions of the Dirac-Brueckner approach and the phenomenological interac- 
tions. The theoretical predictions are also compared with the empirical parametrizations of neutron 
and proton optical-model potentials based on the experimental nucleon-nucleus scattering and the 
phenomenological ones adopted in transport-model simulations of heavy-ion collisions. The direct 
contribution of the rearrangement term due to three-body forces to the single particle potential and 
symmetry potential is discussed. 

PACS numbers: 26.60.-|-c,21.30.Fe,21.65.+f 
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I. INTRODUCTION 



The study on the role of isospin degree of freedom 
is in rapid progress in both nuclear physics and nu- 
clear astrophysics. The experimental and theoretical 
research on isospin physics have received a strong boost 
because of the construction of more advanced detec- 
tors (such as Magnex[3) and new radioactive ion beam 
facilities (project RIA(2|). A wide range of rich phe- 
nomenologies from nuclei far from the /3-stability line 
to strongly asymmetric compound systems formed in 
heavy ion collisions (HIC) requires a deep understand- 
ing of the isospin dependence of the in-medium nuclear 
effective interaction in a large range of nucleon den- 
sity and energy. Among the interesting new physics, a 
key point is the interplay between the isospin T — 
and T = 1 components of the effective interaction as a 
function of the isospin asymmetry. 

Because of the lack of enough empirical informa- 
tion, the most reliable theoretical tools are the mi- 
croscopic parameter-free approaches based on realistic 
nucleon-nucleon (NN) forces fitting the experimental 
phase shifts of the in-vacuum nucleon-nucleon scatter- 
ing. One of the most advanced approaches is the non- 
relativistic Brueckner theory. Over the last decade, in 
fact, it has been improved in two aspects: its con- 
vergence has been verified at the level of three-body 
correlations'^ and the empirical saturation point has 
been reproduced by including microscopic three-body 
forces (TBF) U . Extending the Brueckner calculations 
to spin and isospin asymmetric nuclear matter impor- 



tant predictions have been made on physical quantities, 
including the symmetry energy and the spin and spin- 
isospin Landau parameters[3y . 

Intimately related to the effective interaction is the 
nucleon self-energy, which brings important informa- 
tion on the momentum dependence of nuclear mean 
field, effective mass and optical potential. In the 
Brueckner-Hartree-Fock (BHF) approximation the self- 
energy takes into account not only the interaction of a 
nucleon with inert core (pure BHF mean field) but also 
with core excitations 0, The latter is crucial for 
the comparison with the experimental nuclear poten- 
tial and the optical model potential at low energy. 

The isospin splitting of the nucleon self-energy is the 
main topic of the present investigation. It has been 
calculated in wide ranges of isospin asymmetry, density 
and energy for the sake of application in transport simu- 
lations of HIC as well as for structure calculations close 
to the neutron and proton drip lines. Despite some 
results existing in the literature since long time, the 
present work has been stimulated not only by the new 
opportunity that we can now study the effects of three- 
body forces in the Brueckner theory, but also by the ap- 
pearance of relativistic Dirac-Brueckner-Hartree-Fock 
(DBHF) calculations H [l^. Consistent microscopic 
predictions could provide strong constraints for phe- 
nomenological approaches, which are still affected by 
large uncertainties. These constraints may lead to the 
need of new parametrizations of the Skyrme-like inter- 
actions for the calculations far from the beta-stability 
line. 
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The isovector part of the neutron and proton sin- 
gle particle (s.p.) potentials, i.e., the symmetry poten- 
tial, is one of the basic inputs of the transport models 
for the collisions of radioactive nuclei. In general, the 
shape of the symmetry energy as a function of density 
is determined simultaneously by both the momentum 
and density dependence of the symmetry potential 11], 
therefore the determination of the momentum depen- 
dence of the symmetry potential is crucial for constrain- 
ing the high density behavior of symmetry energy. In 
the earlier dynamical simulations of HIC, the momen- 
tum dependence of the symmetry potential was seldom 
taken into account. Only recently. Das et al. |0| has 
proposed some simple phenomenological parametriza- 
tions for the momentum dependence of the symmetry 
potential which have been ado pted in the dynamical 
simulations of HIC by Li et aL [H El [H where it is 
shown that the experimental observables such as the 
neutron-proton differential flow, the isospin fractiona- 
tion and the 7r~/7r"'" ratio etc., are quite sensitive to 
the momentum dependence of the symmetry potential. 
Microscopically, the proton and neutron s.p. poten- 
tials and their isospin dependence have been studied 
in Ref. |^, |^ within the BHF approach. However, the 
momentum and density dependence of the symmetry 
potential was not discussed in our previous investiga- 
tions jEQ. 

In the present paper, we will concentrate on the dis- 
cussion of the isospin splitting of the effective mass 
and the density and momentum dependence of the 
symmetry potential, based on the BHF approxima- 
tion 0. Especially we shall compare our microscopic 
symmetry potential with the phenomenological ones of 
Ref. T^l and the predictions of the Dirac-Brueckner 
method 9, 10] . The present paper is arranged as fol- 
lows. In Sec. II we present a brief introduction of 
the Brueckner-Bethe-Goldstone (BBG) theory for G- 
matrix with a microscopic TBF, including self-energy 
and effective mass of protons and neutrons. The numer- 
ical results for the symmetry potential are reported and 
discussed in Sec. HI in comparison with other model 
predictions. Summary and conclusions are drawn in 
Sec. IV. 



II. SELF-CONSISTENT BHF APPROACH 
INCLUDING A THREE-BODY FORCE 



A. BBG equation 



The Brueckner theory and its extension to include 
TBFs are described elsewhere 0| . Here we simply give a 
brief review for completeness. The starting point of the 
BHF approach is the reaction G-matrix, which satisfies 
the following isospin dependent Bethe-Goldstone (BG) 



equation, 
G(p,/3,c^) 



(1) 



vnn 



lfcifc2)Q(fcl,fc2)(fclfc2l 



where ki = (ki,ai,Ti), denotes the single particle mo- 
mentum, the z-component of spin and isospin, respec- 
tively, vnn is the realistic nucleon-nucleon (NN) in- 
teraction, Lo is the starting energy. The asymmetry pa- 
rameter is defined as P — {pn — Pp)/p, where p, Pn, and 
Pp denote the total, neutron and proton number den- 
sities, respectively. For the NN interaction, we adopt 
the Argonne Vis two-body interaction ^1] plus a mi- 
croscopic three-body force (TBF) [13. The TBF is 
constructed by using the meson-exchange current ap- 
proach E3 and the most important mesons, i.e., tt, 
p, a and ui have been considered p^ . The parame- 
ters of the TBF model have been self-consistently de- 
termined so to reproduce the Argonne Vig two-body 
force using the one-boson-exchange potential model |j] . 
Their values can be found in Ref. |j]. The TBF con- 
tains the contributions from different intermediate vir- 
tual processes such as virtual nucleon-antinucleon pair 
excitations, and nucleon resonances ( for details, see 
Ref. E3)- The TBF effects on the equation of state 
(EOS) of nuclear matter and its connection to the rel- 
ativistic effects in the DBHF approach have been re- 
ported in Ref. *?] . 

In solving the BG equation for the G-matrix, the con- 
tinuous choice for the auxiliary potential is adopted 
since it provides a much faster convergence of the hole- 
line expansion than the gap choice j^]. One advantage 
of the continuous choice is that the auxiliary potential 
has the physical meaning of the mean field felt by a 
nucleon during its propagation between two successive 
scatterings in nuclear medium ^^l- 

The effect of the TBF is included in the self- 
consistent Brueckner procedure along the same lines 
as in Ref. ^3; where an equivalent effective two-body 
interaction v is constructed by weighting the third par- 
ticle in the real TBF by means of the defect function. 
So doing, one avoids the difficulty of solving the full 
three-body problem. The effective two-body interac- 
tion V can expressed in r-space as0 

{rir2\v\r-^r2) = 



df3df3>:(f3')(l - ??(ri3))(l - 77(^^3)) 



X W3(f/f2'f3'lfif2f3)0„(r3)(l - 7?(ri3))(l - 7?(r23)) 

(2) 

where the trace is taken with respect to the spin and 
isospin of the third nucleon. The function ri{r) is the de- 
fect function. Since the defect function is directly deter- 
mined by the solution of the BG eauationp^. it must 
be calculated self-consistently with the G matrix and 
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the s.p. potential J7(fc)[^ at each density and isospin 
asymmetry. It is evident from Eq.© that the effective 
force V rising from the TBF in nuclear medium is den- 
sity dependent. A detailed description and justification 
of the method can be found in Ref . [I7| . 

B. Selfenergy 

In the BHF approximation 8] with the TBF, the 
selfenergy is made of three terms: 

e) = Sb/i/ + Scpoi + (3) 

where the first term is the HF potential with the G- 
matrix as the effective interaction, the second term is 
due to the core polarization 0, and the third term 
stems from the density dependence of the effective force 
V, i.e., the TBF rearrangement term v |2^. The first 
two terms in the BHF approximation in asymmetric nu- 
clear matter have been discussed elsewhere . In gen- 
eral, the TBF effect on the selfenergy within the BHF 
framework is twofold. First, it affects the selfenergy via 
the modification of the G-matrix. This effect has been 
embodied in the BHF selfenergy, i.e., the first two terms 
in Eq.©. Second, the density dependence of v will in- 
duce an additional contribution, i.e., a rearrangement 
contribution to the selfenergy (third term in Eq.(3)). 
The main contribution of the TBF rearrangement can 
be formally written in the BHF approximation 

ij 

where rii is the Fermi step function. For most results we 
present below the selfenergy is calculated in the BHF 
approximation, and the effect of the TBF is only re- 
stricted to the G-matrix via the BBG equation with 
the two body force and the effective TBF. At the end 
of Sec. Ill we will discuss explicit effect of the rear- 
rangement term of the effective three-body force on the 
self-energy and the symmetry potential. 

When calculated on the energy shell the self-energy 
gives rise to the single nucleon potential. The contribu- 
tion due to the inert core (BHF) for neutrons and pro- 
tons is reported in Fig^ as a function of momentum 
k for three densities and several isospin asymmetries 
f3 = {pn — Pp)/ P- The core polarization term mainly 
influences the potential at k < fc^H ^-nd if "^iH be ne- 
glected in the discussion of this subsection. In order 
to explore the isospin effects on the nucleon effective 
masses (Sec II. C) we split the neutron and proton s.p. 
potentials into the contributions from the isospin T — 
and T = 1 channels, i.e., 

C/p(fc, 13) = Up{k, i3)t=o + Up{k, /?)r=i (5) 
f/„(fc,/3) = [/„(fc,/3)T=o + C/„(fc,/3)T=i (6) 

The isospin behavior of the neutron or proton s.p. po- 
tential is a result of the competition between the T = 



and T = I isospin channels. As discussed in Ref. [3.l2l| 
the isospin effect on the EOS of asymmetric nuclear 
matter is dominated by the isospin T = component 
of the NN interaction. In Fig. [3 we present the contri- 
butions from the isospin T = and T = 1 channels to 
the proton and neutron s.p. potentials at A; = 0, sepa- 
rately, as a function of P with respect to their values in 
symmetric nuclear matter {[3 — Q). It is seen that the 
variations versus P of the T = components are much 
larger than the corresponding T ~ 1 components, i.e., 
3 ~ 5 times larger, implying that the /? dependence of 
the neutron and proton s.p. potentials is determined 
to a large extent by the T = component. This is 
what expected since as increasing the neutron excess, 
the T — interaction between two unlike-nucleons ( 
vanishing between two like-nucleons) becomes stronger 
for protons and weaker for neutrons. The relatively 
small deviations of the T — 1 components of the t/„ 
and Up from their common values in symmetric matter 
is associated to the variations of the Fermi surfaces in 
neutron-rich matter. It is also seen from the figure that 
the net contribution of the T = channels stems al- 
most completely from the SD tensor channel (squares) 
which is strongly attractive at relatively low energies 
while the contributions from other T = channels can- 
cel out each other. This is in agreement with the previ- 
ous observation for nuclear symmetry energy |2lJ |. As 
a consequence, at low momenta, the proton s.p. poten- 
tial becomes more attractive and the neutron one more 
repulsive going from symmetry nuclear matter (/? = 0) 
to pure neutron matter (/3 = 1), as shown in Fig. ^ 
According to the experimental data on the phase shifts 
of nucleon-nucleon scattering, the attraction of the SD 
channel decreases with energy, so that for a given en- 
ergy the T = channel contribution to the splitting 
Un — Up becomes equal to the T — I channel contribu- 
tion. As a result the isospin sphtting {Ut=o — Ut=i)/P 
vanishes for a given value of momentum as marked by 
the crossing point in Fig. ^ This point is almost in- 
dependent of both isospin and density. Therefore, the 
increase of the proton potential depth vs. asymme- 
try results in an increase of the slope as a function 
of momentum. This behavior controls the proton and 
neutron effective mass splitting in neutron rich matter. 
From the above discussion it may be concluded that 
the isospin behavior of the momentum dependence of 
the proton and neutron s.p. potentials which deter- 
mines the neutron-proton effective mass splitting, is es- 
sentially controlled by the tensor component of the NN 
interaction, or say, by the nature of the NN interaction. 
The core polarization affects mainly the s.p. potential 
in the low momentum range below the Fermi surface. It 
gives a repulsive contribution for both the proton and 
neutron potentials and weakens the momentum depen- 
dence at low momentum 0, Q • The core polarization 
contributions to the proton and neutron potential may 
cancel out with each other and thus it modifies only 
slightly the symmetry potential at low momentum 
Therefore inclusion of the core polarization will not al- 
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ter our above discussion and conclusion. 



FIG. 1: Neutron and proton BHF mean fields at different 
isospin asymmetries for three different nucleon densities 



FIG. 2: Isospin splitting the proton and neutron s.p. poten- 
tials at fc = 0; isospin T = (solid curves) channels, T — 1 
(dashed curves) channels, SD tensor channel (squares). 



C. Effective mass 

The nucleon effective mass m* stems from the non- 
local nature of the s.p. potential felt by a nucleon 
propagating in nuclear medium. It is determined by 
the slope of the real part of the on-shell self-energy 
U{k,ek) = ReS(fc,efc) in momentum space, i.e. 



m 
m 



1 



dL/(fc,efc) 
de. ■ 



(7) 



The effective mass itself is momentum dependent, but 
usually its value at the Fermi momentum is considered 
(hereafter we only discuss the latter). The s.p. en- 
ergy is determined by the following momentum-energy 
relation, 



fc2 

Efe = h ReS(fc, Efe) 

2 m 



(8) 



for a given approximation of the self-energy. It is clear 
that the effective mass arises from both the momen- 
tum and energy dependence of the microscopic s.p. po- 
tential. It is a different physical quantity from the 
Dirac effective mass associated to the medium mod- 
ification of the Dirac spinor, which in fact is due 
to the self-consistent requirement in the relativistic 
Dirac-Brueckner approach as it has been clarified in 
Refs. 0, and discussed afterwards. 

In isospin asymmetric nuclear matter since the 
momentum-dependence of the neutron s.p. potential is 
different from that of the proton one, the common value 
of the neutron and proton effective mass is expected to 
split into two branches as a function of isospin asym- 
metry parameter. The BHF result with the TBF con- 
tribution is displayed in Fig. 13 Two main features are 

FIG. 3: Isospin splitting of the effective mass from a Brueck- 
ner calculation with three-body force. 

displayed by the results: the linear dependence on the 
asymmetry parameter /? extended to the whole asym- 
metry range, which is a well-known effect |2ll|. and the 



isospin splitting with m* > m* in neutron-rich mat- 
ter. The latter property is related to the increasing 
(decreasing) slope of the proton (neutron) potential as 
discussed in the preceding subsection, but it has not 
yet received a general consensus. In microscopic theo- 
ries such as the BHF method and the DBHF approach, 
the predicted neutron-proton effective mass splitting 
is such that m* > m* at Fermi surface in neutron- 
rich matter. However, within phenomenological models 
some parametrizations of the extended Skryme-like ef- 
fective interaction may lead to an opposite isospin split- 
ting 0, |23| • A complete settlement of such a contro- 
versy will be one of the most important goals of isospin 
observables such as the neutron-proton differential ffow 
and the ■k~ /ir'^ ratio in HIC induced by radioactive 
beams. 

In order to get a deeper theoretical insight into the 
neutron-to-proton splitting of the effective mass, we 
want to see how it comes out within the Brueckner 
many- body theory. 



FIG. 4: Neutron and proton mass vs. asymmetry: separate 
contributions from the E-mass and k-mass in units of the 
bare nucleon mass. The baryonic density is p = 0.17/m~^ 



The off-shell values of the self-energy I](fc, e) depend 
separately on the energy and momentum and, as a con- 
sequence, following Mahaux et al. 0, the effective mass 
can be split into a product of the two contributions: 



m nif, 
m mm 



where 



me{k) = m 



mk{k) = m 



1 - 



c)S](fc,£) 

de 



1 + 



m dY.{k, e) 



dk 



(9) 



(10) 



(11) 



The fc-mass mk is related to the nonlocality of the mi- 
croscopic mean field in r-space. If the self-energy is 
energy independent (static limit), then m^ = m and 
the fc-mass is equal to the effective mass m*. The e- 
mass describes the nonlocality in time and is related to 
the quasi-particle strength by me{k) = m/Z{k), which 
gives the discontinuity of the momentum distribution 
at the Fermi surface, and measures the amount of corre- 
lations included in the considered approximation. The 
two components of the effective mass for protons and 
neutrons are plotted in Fig.^] It is seen that the isospin 
splitting is opposite for the two effective masses mk 
and me- In neutron-rich matter the neutron fc-mass 
becomes larger than the proton one, while the neutron- 
proton e-mass splitting is reversed. Although the ab- 
solute strengths of the splitting for the two masses are 
about the same, the relative splitting of the fc-mass is 
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much more pronounced than that of the e-mass due to 
the smaher value of the fc-mass. As a consequence, the 
fc-mass predominates the isospin splitting of the total 
effective mass m* and leads to the result m* > m* 
in neutron-rich nuclear matter. This result indicates 
that the effective mass splitting is dominated by the 
nonlocality of the microscopic s.p. potentials in spatial 
space. 

The effective mass so far discussed is also named 
Schrodinger mass to distinguish from the Dirac mass 
appearing in the relativistic mean field theory (RMT) 
and the relativistic DBHF approach ,2^ ^] . The 
Dirac effective mass stems from the self-consistency 
requirement between the s.p. wave function and the 
s.p. spectrum of the Dirac spinor which is dressed in 
nuclear medium, and has no any counterpart in the 
non-relativistic limit. As discussed in Ref. [2^, the 
Dirac mass can be traced back to the effect of the vir- 
tual nucleon-antinucleon pair excitations, i.e., the in- 
medium positive-energy spinor is an admixture of the 
free negative- and positive-energy spinors. In Ref. j22l | , 
the origin of different effective masses defined in liter- 
atures has been discussed in more details. It is shown 
that in the relativistic framework, an nonrelativistic- 
type of effective mass can be introduced based on the 
corresponding Schrodinger equivalent s.p. potential 
and it can be compared to the empirical value ex- 
tracted from analyses in the framework of the nonrel- 
ativistic optical and shell models. Investigations for 
determining the relativistic effective mass in terms of 
the momentum dependence of the s.p. energies in the 
DBHF framework have been made for symmetric nu- 
clear matter in Ref. 12^. Very recently. Van Dalen 
et al. [2^ have investigated the isospin splitting of 
the nonrelativistic-type effective mass obtained from 
the Schrodinger equivalent s.p. potential within the 
DBHF approach. They find that both the dynamical 
structure (i.e., the momentum and density dependence) 
and the splitting of the nonrelativistic-type effective 
mass are satisfactorily consistent with the predictions 
of the nonrelativistic BHF approach. The Dirac mass 
in neutron-rich matter shows an opposite isospin split- 
ting of ^ < mjj p, indicating that the virtual pair 
effect on a neutron spinor becomes stronger as the mat- 
ter goes to neutron-richer. The nonrelativistic-type of 
effective mass in neutron-rich matter derived from the 
RMT displays the same behavior of isospin splitting as 
the Dirac mass, which is due to the fact that the nonlo- 
cal structure of the self-energy is neglected in the RMT 
as discussed in Ref. [g^. 



III. SYMMETRY POTENTIAL AND 
OPTICAL-MODEL POTENTIAL 

The microscopic self-energy in the BHF approach is 
nonlocal in space-time coordinates and thus depends 
on both momentum and energy. When evaluated on 
the energy shell, it corresponds to the empirical optical 



model potential l23| : 

UoptiE)^^ik{E),E), (12) 

where E is the incident energy, and T^{k{E), E) the 
on-shell self-energy. The momentum fc = k{E) is deter- 
mined by the mass-shell relation 

E=^^^ + nk{E),E). (13) 



FIG. 5: BHF isospin symmetry potential vs. momentum k 
for three values of density. 

The isovector part of the s.p. potential, which drives 
the isospin splitting of the nucleon mean field in asym- 
metric nuclear matter, is linearly decreasing with f3 and 
thus the symmetry potential can be defined as 

r T Un Up / -, . \ 

Ufiyra = (14) 

where C/„ and Up are the s.p. potentials felt by a neu- 
tron and proton in nuclear medium, respectively. In 
Fig. O the BHF symmetry potential Ugym is displayed 
as a function of momentum fc for three densities and 
several isospin asymmetries. The nucleon- nucleus scat- 
tering is not influenced so much by the core polarization 
and it can be neglected. It is seen that the symme- 
try potential depends strongly upon both density and 
momentum. Above the Fermi surface [/gym decreases 
rapidly as a function of momentum and saturates at 
high enough momenta. In the momentum region rel- 
evant to the intermediate HIC up to a beam energy 
about 300 MeV per nucleon J7sym is positive, implying 
that its effect is repulsive on neutrons and attractive 
on protons. At higher densities the repulsion (attrac- 
tion) on neutrons (protons) becomes stronger. In Fig. 5 
it is also shown that C/sym is almost independent of the 
isospin asymmetry /3, which indicates that the linear 
dependence of the neutron and proton s.p. potentials 
on [3 persists at any energy and thus it provides a mi- 
croscopic support of the empirical assumption of the 
Lane potential [2^ . 

Experimentally the strength of the Lane potential 
and its momentum dependence can be extracted from 
the nucleon-nucleus scattering data and/or (p,n) charge 
exchange reactions. Earlier optical model analyses of 
the experimental data with incident energies between 7 
and 100 MeV indicate that the {Un — Up)/2(3 at normal 
nuclear matter density has a value of about 28 ± 6 MeV 
at fc = and decreases as a function of incident energy 
with a slope between 0.1 to 0.2pll29l|. 

The Lane potential is represented by the dashed 
results extracted from the experimental data of 
nucleon-nucleus scattering based on the optical poten- 
tial model [25. This area is crossed by both the BHF 
symmetry potential for /3 = 0.2 and the bulk contribu- 
tion of the empirical one from Ref. 30] discussed later. 
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The predicted strength of the isospin sphtting aX k = 
is about 25 MeV in good agreement with the empirical 
value 22 ^ 34 MeV, extracted from the experimental 
data of nuclcon-nucleus scattering based on the optical 
potential model [s^l • 

FIG. 6: Comparison among different symmetry potentials: 
the dashed area is the Lane potential a — bEtin with a — 22— 
34 MeV and b = 0.1-0.2(see Ref. J3l] for more details), the 
dashed line is the BHF result for asymmetric nuclear matter 
with p = 0.16 fm~^ and /3 — 0.2 and the solid curve is from 
empirical parametrization of nucleon-^''*Pb scattering|30|| . 

A recent investigation js^l on a broad range of mass 
(24 < ^ < 209) and incident energy (1 KeV < E < 
200 MeV) provides a new parametrization of the opti- 
cal model potential in terms of volume, surface, spin- 
orbit and Coulomb contributions. Therefore, it should 
be more suitable than the one of Fig. |H| for a com- 
parison with the microscopic potential of nuclear mat- 
ter, including the density, isospin and energy depen- 
dence. Since a direct comparison does not take into 
account the density variation in the nuclear surface, we 
have folded the nuclear-matter mean field with the den- 
sity profile of different nuclei obtained from a Thomas- 
Fermi approximation '2^. So doing we can compare 
the folded mean field with the optical potential as a 
function of /3 = {N — Z)/A. The results are plotted 
in Fig. [71 for incident energy equal to the Fermi en- 
ergy. The theoretical predictions slightly overestimate 
the empirical ones, but the isospin shift turns out to be 
nicely reproduced. Actually the Thomas-Fermi approx- 
imation works better with heavy nuclei for which the 
comparison looks much better. The momentum depen- 
dent optical-model potential of Ref. 30] gets out of the 
dashed area at low energy, indicating that some uncer- 
tainties still affect the optical-model parametrizations. 



FIG. 7: Isospin shift of empirical optical potential at the 
Fermi energy (volume term only) for several nuclei(symbols) 
in comparison with the theoretical nuclear-matter BHF pre- 
dictions (big triangles joined by dashed lines). Small sym- 
bols are the Optical potential fit by Koning and Delaroche 
(NPA 2003). 



FIG. 8: Comparison of the BHF symmetry potential with 
other predictions (see text). 



At high density and energy one may have to rely 
on the scarce information from HIC induced by high 
energy radioactive beams. In this regard, it is interest- 
ing to mention that isospin diffusion has been found 



to be rather sensitive to the momentum-dependence 
of Usym 14]. Up to now only the phenomenological 
parametrizations of the momentum-dependent symme- 
try potential has been adopted in the dynamical simu- 
lations of heavy ion collisions, therefore it is instruc- 
tive to make a comparison between the present mi- 
croscopic symmetry potential with the phenomenolog- 
ical onesjij. In Fig. |S1 is plotted the symmetry po- 
tentials versus momentum for three values of density 
p = 0.085,0.17 and 0.34 fm^''. In the figure the curves 
with filled symbols are the results from the BHF cal- 
culations, repeated in the four panels, while the curves 
with open symbols are the phenomenological ones of 
Ref. ^31 and different panels correspond to different 
parametrizations. It is clear from the figure that the 
microscopic J7sym shows a remarkably different behavior 
from the phenomenological ones as a function of den- 
sity and momentum. All of the four phenomenological 
symmetry potentials drop much faster at high density 
p = 0.34 fm""^ as the momentum increases as com- 
pared to our BHF one. In the two cases of the GBD(O) 
and GBD(l) parametrization, the deviation from the 
BHF prediction is especially large in the whole density 
and momentum regions considered here. At the nor- 
mal nuclear density, the momentum dependence of the 
Gogny and MDI(O) parametrization is closer to our mi- 
croscopic one, but the Gogny Ugym presents an opposite 
density dependence and discrepant dramatically with 
the microscopic one at high densities. For example, at 
p — 0.34 fm~'^ the Gogny Usym is attractive (repulsive) 
for neutrons (protons), while the BHF one is repulsive 
(attractive) for neutrons (protons) up to fc ~ 4 fm~^. 
Even the MDI(O) parametrization which is closest to 
our microscopic prediction, turns out to become quite 
different at high densities. 

As the last point we discuss the effect of the rear- 
rangement contribution of the TBF. Due to the density 
dependence of the effective force v, the TBF provides an 
extra repulsive contribution Stb/ (see Eq.®) to both 
the proton and neutron s.p. potentials. At high density 
this contribution is expected to be strongly momen- 
tum dependent and may affect considerably the high 
momentum components of the fragmentation residues 
in HIC. We compare our calculated results with the 
parametrization of the optical potential used in the 
trans por t-model simulations of elliptic flows in central 
HIC I33, where high densities are reached. The po- 
tentials in symmetric nuclear matter at p = 0.3fm~'^ 
are shown in Fig El where the line with full squares 
is the one of Ref. ^32] which has been shown to de- 
scribe the observed elliptic flow data fairly well. It is 
clear from the figure that the BHF potential without 
the TBF is too attractive, especially at hi gh d ensities, 
as compared to the one proposed in Ref. p3^. and its 
momentum dependence at high momenta turns out to 
be too weak [33 for describing the experimental ellip- 
tic flow data. Inclusion of the TBF effect only via the 
G-matrix, i.e., inclusion of the TBF effect in the first 
two terms of Eq.©, weakens the dependence of the 
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s.p potential on momentum,^]. It is seen from Fig. |3 
that the rearrangement contribution of the TBF, i.e., 
the third term Tm, f of Eq. Q , provides a strongly extra 
repulsion to the optical potential and improves remark- 
ably the agreement between our microscopic potential 
and the parametrized one of Ref. [s^l ■ We find that the 
TBF rearrangement leads to also a strong momentum 
dependence at high densities. For instance at p = 0.3 
fm^'^ the extra term turns out to be AU = 5.68-|-6.84A;^ 
in units of MeV. 

To study the effect of the TBF rearrangement on the 
isospin symmetry potential, we report in Fig. llOl the re- 
sults in coniparison with the recent DBHF predictions 
from Refs. Therein it is displayed the neutron 

s.p. potential (left-upper panel), the proton s.p. po- 
tential (left-lower panel) and the symmetry potential 
(right panel) in neutron-rich matter with P = 0.4 for 
both cases with and without including the Titb f term. It 
is seen that the BHF neutron and proton s.p. potentials 
without the ^tbf term are much more attractive than 
the DBHF ones. Inclusion of the Titbf term leads to a 
strong enhancement of the repulsion of both the proton 
and neutron s.p. potentials and reduces substantially 
the disagreement between the s.p. potentials predicted 
by the nonrelativistic BHF and the relativistic DBHF 
approaches. It is shown in the right panel of Fig. ^| 
that the effect of the Htbf term on the isospin sym- 
metry potential is very weak, indicting that the con- 
tributions of the Sth/ term to the neutron and proton 
potentials cancel out in a wide momentum range with 
each other to a large extent. In both cases with and 
without the Stb/ term, the isospin symmetry poten- 
tials obtained by the BHF approach display an over- 
all agreement with those by the DBHF approach. At 
high momentum our symmetry potential is slight lower 
than the DBHF ones. The difference between the two 
DBHF calculations may be attributed to the different 
methods adopted to extract the selfenergy as discussed 
in Ref. where it is shown that the determination of 
the nucleon selfenergy in the DBHF framework is still 
affected by some uncertainties. 



FIG. 9: Optical potential in nuclear matter as a function 
of nucleon energy at density p = 0.3/m~^ from the BHF 
calculation including all the three terms of EqElwith the 
Argonne Vis interaction plus TBF (solid line). The squares 
are from Ref. . The circles are from the earlier BHF 
calculation without three-body force [s^ . 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have investigated the momentum 
and density dependence of the symmetry potential and 
discussed the origin of the neutron-proton effective 



FIG. 10: Comparison of the BHF symmetry potential with 
and without the T^ttf contribution with Dirac-Brueckner 
predictions. The curve with full circles is from Ref. H^ . 
the one with full squares is from Ref.Q 



mass splitting in neutron-rich nuclear matter within the 
framework of the Brueckner theory. 

We have found that the isospin behavior of the mo- 
mentum dependent neutron and proton s.p. potentials 
can be traced back to the effect of the SD tensor com- 
ponent of the NN interaction and, consequently, the 
neutron-proton effective mass splitting is essentially de- 
termined by the intrinsic properties of the NN interac- 
tion. The obtained neutron-proton effective mass split- 
ting in neutron-rich matter is m* > m* in good agree- 
ment with the recent predictions by the non relativis- 
tic limit of DBHF approach . The isospin splitting 
TO* > m* is shown to stem from the splitting of the 
/c-mass, i.e., from the spatial nonlocality of the micro- 
scopic neutron and proton s.p. potentials. 

It turns out that the predicted symmetry potential 
depends sensitively on density and momentum, but 
almost independent of the isospin asymmetry (Lane 
ansatz). In the energy and density regions most rele- 
vant for the nucleus-nucleus scattering up to an incident 
energy of about 300 MeV per nucleon, our microscopic 
symmetry potential is repulsive for neutrons and at- 
tractive for protons, and its strength becomes smaller 
as momentum increases for a fixed density. A satisfac- 
tory support to the microscopic predictions is provided 
by a recent optical model parametrizations of nucleon- 
nucleus scattering [s^l- In dense nuclear matter, which 
can be probed in HIC at intermediate and high energies, 
the symmetry potential turns out to become stronger 
in the high momentum region up to about 4 fm~^ as 
increasing density. At high energy relativistic effects 
manifest with a strong momentum dependence |26ll32j |. 
which can be interpreted as an effect of the TBF. At 
the present time a comparison with phenomenological 
Skyrme-like or Gogny predictions of the symmetry po- 
tential |3lj is quite difficult, and empirical constraints 
on their parameters are needed to make it useful. 

In the BHF approach, the TBF contribution has been 
included by reducing the TBF to an equivalent effec- 
tive two-body force. The rearrangement effect due to 
the density dependence of the equivalent force has been 
found to provide an extra repulsive contribution to the 
proton and neutron s.p. potentials, which improves 
substantially the agreement of our nonrelativistic s.p. 
potential with the parametrized potential for describ- 
ing the elliptic flow data|33l and those predicted by 
the DBHF approach llO^. The TBF rearrangement 
has been found to affect only slightly affect the isospin 
symmetry potential due to the cancellation between the 
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two contributions to the proton and neutron potentials. 
Our calculated symmetry potential is shown to be in an 
overall agreement with the DBHF predictions. 

Physical observables which are sensitive to the sym- 
metry potential, including the neutron-to-proton ratio 
of pre-equilibrium nucleon emission, neutron-to-proton 
differential flow and isospin diffusion, are expected to 
provide experimental constraints on the momentum 
and density dependence of the symmetry potential |3lj | . 
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